Statistical properties of zeta functions' zeros 



o 

(N 



(N 



V. Kargin 

282 Mosher Way, Palo Alto, CA 94304 
e-mail: Vladislav. karginSgniail . com 



Abstract: The paper reviews existing results about the statistical distribution of zeros 
for three main types of zeta funetions: number-theoretical, geometrical, and dynamical. 
The paper provides necessary background and some details about the proofs of the main 
' results. 

D 

I ' AMS 2000 subject classifications: 11M26, 11M50, 62E20. 

Keywords and phrases: Riemann's zeta, Selberg's zeta, Ruelle's zeta, Montgomery's 
conjecture, distribution of zeros. 



PLh , 1- Introduction 

^- 

. The distribution of zeros of Riemann's zeta function is one of the central problems in modern 

mathematics. The most famous conjecture is that all of the zeros are on the critical line 
Res ~ 1/2. A more precise version says that the zeros behave like eigenvalues of large Hermitian 
matrices. While these statements are conjectural, a great deal is known about the statistical 

■ properties of these zeros and about zeros of closely related functions. In this report we aim to 
summarize findings in this research area. 

I We give necessary background information, and we cover the three main types of zeta 

functions: number-theoretical, Selberg-type, and dynamical zeta functions. 
' Some interesting and important topics are left outside of the scope of this report. For 

■ example, we do not discuss quantum arithmetic chaos or characteristic polynomials of random 
matrices. 

The paper is divided in three main sections according to the type of the zeta function we 
' discuss. Inside each section we tried to separate the discussion of the properties of zeros at the 

H ' global and local scales. 

■ ■ ■ Let us briefly describe these types of zeta functions and their relationships. First, the 

number-theoretical zeta functions come from integers in number fields and their generaliza- 
tions. Due to the additive and multiplicative structures of the integers, and in particular due 
to the unique decomposition in prime factors, the zeta functions have a functional equation, 
C (1 - s) = c (s) C (s) and the Euler product formula C (s) = Dp (1 - (^P)")"^ • 

One can also look at number-theoretical zeta functions from a different prospective if one 
starts with modular forms, that is, functions on two dimensional lattices that are invariant 
relative to transformations from SL2 (Z) . 

Some modular forms (the cusp forms) can be written as X]n>o '^■'^P (27rmz) , where z 
is the ratio of the periods of the lattice, and one can associate a zeta function ^ c„n~''' to 
this modular form. Then the modularity ensures that the zeta function satisfies a functional 
equation. In addition, the zeta function will have the Euler product property if the original 
modular form is an eigenvector for certain operators (the Hecke operators). 
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This class of modular zeta functions overlaps with number-theoretical functions to a large 
extent. 

In fact, a recent significant development occurred when it was proved that all zeta functions 
associated with elliptic curves come from modular zeta functions. Among other applications, 
this discovery was a key to the proof of Fermat's last theorem. 

An important representative of the second class of zetas is Selberg's zeta function. Let H 
is the upper half-plane with the hyperbolic metric and F is a discrete subgroup of SL2 (Z). 
Selberg showed that sums over eigenvalues of the Laplace operator on Riemann surface F\H 
are related to sums over non-conjugate elements of F. This relation is called Selberg's trace 
formula. Selberg's zeta function is constructed in such a way that it has the same relation to 
the trace formula as Riemann's zeta function has to the so-called "explicit formula" for sums 
over zeros of the zeta function. 

While Selberg's zeta function resembles Riemann's zeta in some features, there are signifi- 
cant differences. In particular, the statistical behavior of its zeros depends on the group F and 
often it is significantly different from the behavior of Riemann's zeros. 

The third class of zetas, the dynamical zeta functions, appeared first as a generalization of 
Weil's zeta functions. 

Let kq is a finite field with q elements and / {X, Y) is an absolutely irreducible polynomial 
over this field. The simplest Weil's zeta function is the usual number-theoretical zeta functions 
for integer ideals of the field kg [X, Y] / f {X, Y) . It turns out (non-trivially) that this zeta 
function can be written as a generating function for fixed points of powers of the Frobenius 
map: a; — > x"^. (The Frobenius map acts on the points of the curve / {X,Y) in the algebraic 
closure of kg.) 

The theory of dynamical zeta functions starts from this setup and defines a zeta function 
for arbitrary map _F on a set M as the generating function for fixed points of powers of F. 
This definition can also be generalized to fiows on a set, that is, to maps F : M x R+ — > M. In 
particular, Selberg's zeta function can be understood as the dynamical zeta function for the 
geodesic flow on the surface F\H. 

With this overview in mind we now come to a more detailed description of available results 
about the statistics of zeta function zeros. 

2. Number-Theoretical Zetas 
2.1. Riemann's zeta 

There are several excellent sources on Riemann's zeta and Dirichlet L-functions, for example 
the books by Davenport [8] and Titchmarsh [47]. In addition, a very good reference for all 
topics in this report is provided by Iwaniec and Kowalski's book [21]. 
By definition Riemann's zeta function is given by the series 



for Res > 1. It can be analytically continued to a meroniorphic function in the entire complex 




(1) 
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plane and it satisfies the functional equation 



Indeed, we can relate C (s) to the series 9 (x) = J2'^=i e " , 



r(i.)C(-s) 



-s/2 



poo poo 
„-i "'O Jo 



From the identities for the Jacobi theta-functions, implied by the Poisson summation formula, 
it follows that 

2.(.) + l = -L(2.(-L)+l 

Writing 

roo pi poo 

a;"/2-i6' (x) d.T = y x'^'^-^e{x)dx + J x'/'^-^9 (x) dx, 
and applying the identity to the integral from to 1, wc find that 

which is symmetric relative to the change s — > 1 — s. 

The second fundamental property of Riemann's zeta is the Eulcr product formula: 

c(»^n(i-^)" (3) 

valid for Res > 1. 

Let the non-trivial zeros of Riemann's zeta function be denoted by pk = /Jfc + ijk- The 
Riemann's hypothesis asserts that /3a; = 1/2 and it is known that < /3fc < 1. The zeros are 
symmetric relative to the real axis, so it is enough to consider zeros with positive imaginary 
part, 7fc > 0. We order them so that the imaginary part is non-decreasing, 71 < 72 < 

One of the most important ideas in the field is to relate the sums over prime numbers and 
sums over zeta function zeros. For example, the Ricmann-Mangoldt formula says that 

n<x p n 

where A [n) = logp, if n is a prime p or a power of p, and otherwise A [n) = 0. The idea of the 
proof (due to Riemann) is to consider the logarithmic derivative 

^ n=2 



and then to use the following formula: 



ifO<j/<l, 



1 r+'^ sds , , 

1 if y > 1, 
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where c > 0, in order to pick out the terms in a Dirichlet series with n < a; by taking y = x/n. 
From (5), one gets 

, ds 



c+ioo r >/ 



Moving the hne of integration away to infinity on the left and collecting the residues at the 
poles one finds formula (4). (See Chapter 17 in Davenport [8] for a detailed proof.) 
More generally, the Landau's formula holds: 



E 

n<.x 



A{n) 



ry,/-'—^ 



1 — s '■ — ' p — s 

p 



2n + s C 



(6) 



Formula (4) is useful if one want to study the distribution of primes and something is 
known about the distribution of zeros, and the formula (6) can be used to study the behavior 
of ^ (s) if some information about the primes is known. 

One difficulty is that the sums over zeros in these formulas arc only conditionally convergent, 
(pfc w 1/2 + 27rz-j^^). In particular, in these sums it is assumed that the summation is in the 
order of increasing |Imp|. 

Another variant of this formula was discovered by Sclbcrg and it avoids the problem of 
conditional convergence. Define 



in) 



A{n) 



A(n), 

log^ ^-21 



2 log^ X 
ANSI: 



Then, 



is) 



E 



A^ (n) 



log^ X 



+■ 



E 

p 

1 



for 1 < n < X, 
for X < n < x^, 

for x^ < n < x'^. 
xP-" (1 - xP-^f 



{s-pY 



(7) 



log^ X • (1 — s)^ 



log^ X 



E' 



{2q + sY 



Finally, there is one more variant of the explicit formula (sometimes called Delsarte's ex- 
plicit formula). Suppose that H (s) is an analytic function in the strip — c < Im s < 1 + c 



(for c > 0) and that \H {a + it)\ < A{1 + \t\) 
iJ (i + «i) and define 



-(1+5) 



uniformly in a in the strip. Let h {t) 



h (x) = / h{t)e 



-2'iritx 



dt. 



(Note that analyticity of H (s) implies that h (x) has finite support.) Then, 



p 71 — 1 ^ 



logri 
2tt 



log n 
2n 



(8) 



where 



r' /I 



r' /I 
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The idea of the proof is similar to the proof of the Riemann-Mangoldt formula. One starts 
with the formula 



1 



2+200 



(^) 



H (s) ds 



n=l -J -co 

1 A(n)^ /logn\ 

n—l ^ ^ 



Next, one can move the line of integration to (—1 — icx), —1 + ioo) and use the calculus of 
residues and the functional equation to obtain formula (8). 

The most general explicit formula was derived by Weil (see [51]). We will not present it 
here since it involves adelic language and this would take us too far apart. 



2.1.1. Statistics of zeros on global scale 

Let Af{T) denote the number of zeros with the imaginary part strictly between and T. If 
there is a zero with imaginary part equal to T, then we count this zero as 1/2. Define 

5(r) :=ilmlogC(i+ir), 

TT Z 

where the logarithm is calculated by continuous variation along the contour a + iT, with a 
changing from +oo to 1/2. 

By applying the argument principle to C and utilizing the functional equation, it is possible 
to show (see Chapter 15 in [8]) that 

T T 7 /I 
AA(T) = _ log — + ^ + 5 (T) + O ' 



27r °27re 8 " ' \l + T 
Let 

Xit) 

Vloglogl 

Then, we have the following theorem by Selberg. (See [42] and [41].) 

Theorem 2.1 (Selberg) Assume RH, and let T" < H < T^, where a > 0. Then for every 
k > 1 

1 /•^+^, .. 2fc! 

H 

with the constant in the remainder term that depends only on k and a. 



|x(0p''^rfi = ^ + o(i/iogiogr). 



In other words, X [t) behaves like a Gaussian random variable. Note that the Riemann 
Hypothesis is not assumed in this result. 

Recently, this was generalized by Bourgade in [7], who found the correlation structure of 
X{t). 

Theorem 2.2 (Bourgade) Let ur be uniform on (0,1), et 0, 1/logt, and functions 

< /t < . . . < < c < cxj. Suppose that for all i ^ j, 



los 



log et 



-> G [0,00] 
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Then the vector 



^ log C ( J + + + i^t, ...l+£t + ifP + iujt 



\J- log et V V2 ' ■ ' ■ ' 2 



converges in law to a complex Gaussian vector (Yi , . . . , Y; ) with the zero mean and covariance 
function 

\ 1 A Cij if i ^ j. 

Moreover, the above result remains true if et ^ 1/ log t, replacing the normalization — log St 
with log log t. 



Cov(r,,y,: 



This result implies the following Corollary for the zeros of the Riemann's zeta. Let 

-^logi^^ 



A(ii,i2) = (AA(t2)-^log^)-(AA(ti) 



2tt 



27re J ' 



which represents the number of zeros with imaginary part between ti and t2 minus its deter- 
ministic prediction. 

Corollary 2.3 (Bourgade) Let Kt be such that, for some e > and all t, Kt > e. Suppose 
log A' t / log log t — > (5 e [0, 1), as t oo. Then the finite- dimensional distributions of the process 



A {ujt + a/Kt, ujt + P/Kt) 



, < a < /3 < oo, 



converge to those of a centered Gaussian process (A (a, P),Q<a<(3<oo) with the covariance 
structure 

1 if a = a' , and /? = /?', 

1/2 ifa^ a', and /3 ^ /?', 

1/2 ifa^ a', and /3 = /?', 
-1/2 = 
elsewhere. 



E 



(A(a,^)A(a',/3')) = < 



Note that Kt is of order (logi)~ and the average spacing between zeros is 1/logt, hence 
the number of zeros in the interval {ujt + a/Kt, cot + P/Kt) is of order (logt)^ . 

This result perfectly corresponds to the results of Diaconis and Evans in [9] about fluctu- 
ations of eigenvalues of random unitary matrices. 

Selberg and Bourgade's theorems are based on an approximation the functions S{t), dis- 
covered by Selberg. 

Proposition 2.4 (Selberg) Suppose k € Z+, < a < 1. Then there exists Ca.k > such that 
for any 1/2 < cr < 1 and t"-/'' < x < t^^'', it is true that 



logC{a + ts)- J2 i 



2k 



ds < Ca,k 



The proof of this statement, in turn, depends on Selberg's formula (7). 
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2.1.2. Statistics of zeros on local scale 



At the local scale, we have a result due to Montgomery, [32]. Assume the Riemann Hypothesis, 
and define 



° 0<7,7'<T 

where a is real, T > 2, and w (u) is a weighting function, 

^ ' 4 + 1*2 

Theorem 2.5 (Montgomery) The function F (a) is real and F (a) = F {—a) . If T > 
Tq (e) , then F (a) > — e for all a. For fixed a satisfying < a < 1 we have 

F{a) = (l + o(l))T-2"logr + a + o(f), 

as T tends to infinity; this holds uniformly for Q < a < 1 ^ e. 

In fact, this result holds uniformly throughout < a < 1, as later proved by Goldstone. 
This result can be used through the identity 

/ log r \ T roo 

rf(7-7')— — j ■i«(7-7') = — logT / F{a)r{a)da, 

0<7,7'<T ^ ^ •'^^^ 

where 

r{a) = / r (w) e-2^'""du. 



Montgomery conjectures that 

F{a) = l + o(l) 

for a > 1, uniformly in bounded intervals. It can also be formulated as follows. 
Conjecture 2.6 (Montgomery) For fixed a < f3, 



E 



0<7,7'<T 
27rQ/logT<|7-7'|<27r;3/ logT 













/ sin TTu \ ^ 


du 


(/: 








\ TTU J 





T 
2^ 



as T goes to infinity. Here d (a, /3) = 1 i/ G [a, j3] , d (a, /3) = otherwise. 

The proof of Montgomery's theorem is based on the following variant of the explicit formula. 
Lemma 2.7 (Montgomery) If 1 < a < 2 and a; > 1 then 

(logr + (1)) + 0„ (^1/2^-1) , 
where t = \t\ + 2. The implicit constants depend only on a. 
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In a further development Hughes and Rudnick in [17] study the sums / (xi) + . . . f (xn) 
where xjy denote the scaled zeros and calculate their moments. This is possible to do if / is 
sufficiently smooth. They proved the following result. Write the zeros as ^ + i^j where 7^ are 
not assumed real. Let / be a real-valued even function with smooth Fourier transform that 
has compact support. Define 

Mfit,T) ■.^Y.f 

(If / were the indicator function of an interval [—1, 1] and if all 7^ are real, then this function 
would count number of zeros in the interval — 27r/ logT, t + 27r/ logT] . However, the Fourier 
transform of this indicator function does not have compact support.) 

Next, choose a weight function w {x) , such that w > 0, J w (x) dx — 1, and w {x) is 
compactly supported, and define an averaging operator 

Theorem 2.8 (Hughes-Rudnick) Assume that the averaging window H ^ T°- for < a < 

1, and that the support of f {u) ~ j f {x) e~^'^™"(ia; is in (— 2a/m, 2a/ m) with integer m> 1. 
Then, as T ^ 00, the first m moments of J\ff, I^N^^ converge to those of a Gaussian 
random variable with expectation J f (x) dx and variance 

aj^J min {\u\ ,1) f (u)"^ du. 

They note that a similar result holds in the random matrix theory for eigenvalues of a 
unitary random matrix, and that the higher moments converge to non-Gaussian values. This 
is again in agreement with the situation in the random matrix theory. 

We will describe the ideas of their proof below in the case when they are applied to zeros 
of Dirichlet's L-functions. The key idea is to use a variant of the explicit formula (8). 

A generalization of Montgomery results about pair correlations of zeros is due to Rudnick 
and Sarnak in [37] who studied statistical properties of fc-tuples of zeros. They proved these 
results for a quite large class of zeta functions, and we will postpone the discussion of their 
results to a later section. 

Rudnick and Sarnak proved their results under some restrictive conditions, similar to con- 
ditions imposed by Montgomery. Roughly speaking these conditions require that the test func- 
tions used for evaluation of correlation functions are very smooth. It is an outstanding problem 
to prove that all results about correlations of zeros hold without these restrictive hypotheses. 
This is supported by numerical evidence and, if this conjecture true, wc can conclude that 
the behavior of the zeta zeros is very similar to the behavior of eigenvalues of random unitary 
matrices. 

2.2. Dirichlet's L-functions 

Let X C*^) denote a multiplicative character modulo a positive integer q. That is, the function 
X maps integers to the unit circle; it is multiplicative, x (nm) = x ("■) X i'"^) : ^-nd x ('^) = if 
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n and q are not relatively prime. The character which sends every integer relatively prime to q 
to 1 is called the principal character modulo q. The conductor of the character is the minimal 
integer N such that the character is periodic modulo N. For simplicity, let g be a prime in the 
following. In this case the conductor equals q. A character is odd if x (^1) = ^-nd even if 

Dirichlct's L- function corresponding to the character x is defined by the series 



L{s,x) = J2 



n—l 



This series can be continued to a function which is meromorphic in the whole complex plane 
and satisfies a functional equation. Namely, let /.j = if x is even and 1 if x is odd. Define 



(5, x) = qi(^+^)n-^(^+^'>T Q {s + ^i)^ L {s, x) , 



then the functional equation has the form 

$(l-5,x)-^$(s,x), 

where t (x) is the Gauss sum: 

m—l 

Note that |t (x)| = ^/q■ For a proof see Chapter 9 in Davenport 
Finally, note that if x is not principal, then L (s, x) is entire. 



2.2.1. Global scale 

For T > 0, let N (T, x) denote the number of zeros of L (s, x) with < <t < 1 and < t < T, 
possible zeros with t = or t = T counting one half only. Let 

S{t,x) = -lm\ogL{^+it,x)- 

Then it can be shown that 

^' ^ £ ^ - ^ + ^' - ^ <»• ^> + (tTt) ■ 

See formula (1.3) in Selberg's paper [40] and for example. Chapter 16 in Davenport [8]. 

If the character is fixed and the imaginary part of the argument is large, then the results 
for Dirichlct's i-functions arc similar to those for Riemann's zeta function. 

A different situation arises when the interval [0,r] is fixed and the character x is random. 
This situation was studied by Selberg, who proved the following results. (Sec [40].) 

First, for \t\ < g^^"'"^, the difference 

p<q y 

is uniformly bounded in the mean as g — ^ oo, and |S'(f,x)| is of the order of magnitude 
Vlog log g in the mean. Selberg also established the following theorem. 
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Theorem 2.9 (Selberg) For \t\ < q^l^ ^ , we have 

where the summation if over all non-principal characters over the base q. 

In other words the distribution of S {t, x) approaches the distribution of a Gaussian random 
variable with variance log log g. In these results t is assumed to be fixed. 



2.2.2. Local Scale 



Hughes and Rudnick in [18] study the linear statistics of low-lying zeros of L-functions on the 
local scale. Hughes and Rudnick order the zeros pi^^ — i + *7i,x follows: 

. . . < Re7_2,x ^ ^'^l-i-x < — ^-^71^^ < Re72,x < ■ • ■ 

and define 

_ logg 

^i,X ~ 271" 

Then they define 

Wfix)= E /(^^x) 

i— — oo 

where / is a rapidly decaying test function. 

The question is to understand the behavior of the averages 

The basis for their analysis is a variant of the explicit formula (8) that relates a sum over 
zeros of L (s, x) to a sum over prime powers. This formula is a particular version of the formula 
from Rudnick and Sarnak [37], which is valid for zeta functions of all cuspidal automorphic 
representations (which we will not define in this paper). Let 



Ai(x) = 



0, if X is even, 

1, if X is odd, 



and let h (r) be any even analytic function in the strip —c < Imr < 1 + c (for c > 0) such 
that \h (r)| < ^ (1 + |r|)"^^+''^ (for r e A > 0, S > 0). Set h{u)^ h{r)e~'''''dr, so that 
h{r) = Jh{u)e'''"du. Then 

1 f°° 

E^(t^-^-) = ^ h{r){logq + G^{r))dr 

- E ^^'^ ") + X («)) , 



where 



Gx ('') = f ^ + ix) + ir] + [\+ p-ix) -ir] - ^ logvr. 
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By using this formula, Hughes and Rudnick separate Wf (x) into two components. Take 
h{r) ^ f (^^^^^r^ , so that h (u) ~ u^f (logq) • have to impose some conditions on / 
so that the conditions on h (r) arc satisfied. We say that / is admissible, if it is a real, even 
function whose Fourier transform / (u) := J f{r)e~'^'^'^™dr is compactly supported, and such 
that 1/ (r)| < ^ (1 + Irj)"*-^'''''-' . Then, we get the following decomposition: 

where 

Wjix) /^/(^^) i^ogq + G^ (r))dr, 

and 

»'r(..:^-i4i:^/(S?)(.(")-^(..)). 

For large q, 

which is asymptotically independent of x- 

First, Hughes and Rudnick show that if Supp ^/^ C [—2, 2] , then 

Next, they prove the following result. 

Theorem 2.10 (Hughes and Rudnick) Let f be an admissible function and assume that 

Supp (^fj C [-a, a], a> 0. 
If m < 2/a, then the m-th moment of Wj"'^ is 

WmliWr^D =1 2-^^/"w2)! ^(/)"N ^frmseven, 
q^^w s ' 'q y 0, if mis odd, 

where ^ 

aif f= J Julfiufdu. 

They also made a conjecture about the limit when the support of / is not bounded to 
[— 2/to, 2/m] by performing an analogous calculation for unitary random matrices (where all 
moments of the eigenvalue distribution are known from the work of Diaconis and Shahshahani 
[10]). A consequence of this conjecture is that the limit of the linear statistic Wf (x) is not 
Gaussian for high moments. 

Hughes and Rudnick used their results to derive some results for the smallest zero of 
L (s, x) • In particular they show that for infinitely many q there are characters x such that 
the imaginary part of the zero is smaller than 1 /4, and that for all sufficiently large q, a positive 
proportion of L (s, x) has a zero smaller than a certain constant. 
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2.3. L- functions for modular forms (The Hecke L-functions) 

A good introductory source for the material in this section is Chapter V of Milne's book [30] 
about elliptic curves. 

2.3.1. L-functions from modular forms 

Let r is a subgroup of finite index in SL2 (Z) and let H denote the upper half-plane {z|Imz > 0} . 
The set H* = HU {00} U Q can be made into a Hausdorff topological space and one can define 
a continuous action of SL2 (Z) on H* as an extension of the action of SL2 (Z) on H : 

(a b \ . , . . (az + b 
[c d j-^(^) = ^ 

The cusps are points in H*\H. One can also show that r\EI* is a compact Hausdorff space 
(that is, a compact space in which every two points have disjoint open neighbourhoods), which 
is a Riemann surface (that is, it admits a complex structure). Recall the following definition. 

Definition 2.11 A modular form for T of weight 2k is a function H — > C such that 

(a) / is holomorphic on H; 

(b) forany7= ( " |j ) e T, f (jz) = (cz + df f (z) ; 

(c) / is holomorphic at the cusps. 

In particular, analyticity at infinity implies that a modular form can be written as 



/ (z) = ^ c (n) q", where q = e^'^' 



z/h 
n>0 



where h is the minimal integer /i > 0, such that ( I G T. 

y I y 

A cusp form is a modular form which is zero at all the cusps. It can be written 



/(z) = ^c(n)g". 



It is usual to write M^^k (r) and S2k (r) to denote the vector spaces of all modular and 
cusp forms. 

Recall that Fq (TV) is a subgroup of SL2 (Z) defined as follows: 

ro(iV) = |(^ ^ c^0(mod7V)|. 

Definition 2.12 Let f he a cusp form of weight 2k for Tq (N) , 

/(z) = ^c(n)g". 

n>l 

The L -series of the cusp form f is the Dirichlet series 

L{f,s) = Y,c{n)n-^. 
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It is possible to estimate that |c(n)| < Cn , and therefore this series is convergent for 
Res > fc + 1. 

Let aN = I I ■ Then 

\ N J 

-c/N \ 




and so the conjugation by a at preserves Fq (N) . This implies the following fact about cusp 
forms. Define operator wjv by the following formula: 

2k 



Then wjv preserves (To (N)) and wf^ = 1. Hence the only eigenvalues of wn are ±1, and 
(Lq (N)) is a direct sum of the corresponding eigenspaces, S2k = S^/} + 5;^^. 
It turns out that the zeta functions corresponding to these eigenvector cusp forms satisfy 

functional equations analogous to the functional equation of the Riemann's zeta. 

Theorem 2.13 (Hecke) Let f e S2k (To (N)) be a cusp form in the e-eigenspace, e = 1 or 
— 1. Then f extends analytically to a holomorphic function on the whole complex plane, and 
satisfies the functional equation 

AisJ) = ei~l)'A{k~sJ), 

where 

A(s,/) = Ar^/2(2^)-^r(s)L(/,s). 

Moreover, under certain conditions, which we explain below, these zeta functions can be 
written as a Euler product over primes. 

Namely, the Hecke operators T (n) (n > 1) are linear operators on S2k (Fq (N)) . We skip 
their definition here, which can be found in Section V.4 of Milne [30] or Section VII. 5 of Serre 
[43] , and only list their properties. 

Theorem 2.14 (Hecke) The maps T [n) have the following properties: 
(a) T (mn) — T (to) T (n) if gcd (to, n) — 1] 

(h)T (p) T (p'') = T (p''+i) + p2fc-iT (p'-^i) ifp does not divide N; 

(c) T{p^) = T{p)\r>\, ifp\N; 

(d) all T [n) commute. 

Since the Hecke operators commute we can look for the functions which are eigenfunctions 
for all of them. Then, it turns out that these eigenfunctions have Euler products. 

Theorem 2.15 (Hecke) Let f he a cusp form of weight 2k for Tq (N) that is simultaneously 
an eigenvector for all T (n) , say T [n) / = A [n) /, and let 

f = Y,c{n)q\ q = e^^". 

Let c(l) = 1. Then, (i) coefficients of the series are eigenvalues of the Hecke operators, 

c (n) = A (n) 
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and (a) 



- — c(p) p ^ ^ I — c (p) p + n^'"" 1 

In particular, S12 (Fq (1)) has dimension 1, and therefore it is generated by a single function, 
which is called the A-function: 



where t (n) is called the Ramanujan r-function. It follows that the L-function associated to 
the A-function has both a functional equation and the Euler product property. 

In a more general situation, if we wish to find forms that have both a functional equation 
and a Euler product, then we must overcome the obstacle that operators wn and T [n) do not 
always commute. However, this obstacle can be circumvented and it can be proved that such 
modular forms do exist. They are called newforms. 

In summary, the L-functions of modular newforms have both a functional equation and the 
Euler product property. 



2.3.2. L-functions from Maass forms 

A nice source for the material in this section is the lecture notes by Liu [25] . A lot of additional 
information about Maass forms can be found in the book by Iwaniec [20]. 

Modular forms are holomorphic and so are comparably rare animals. It is not immediately 
clear how to compute them. We can try to relax the assumption of complex analyticity, and 
in this way we come to the concept of a Maass form. 

Definition 2.16 A smooth function f ^ is called a Maass form for group F, if 

(i) for all .9 G F and all z G H, / [gz) = f (2) ; 

(ii) f is an eigenfunction of the non-Euclidean Laplace operator: 

and 

(Hi) there exists a positive integer N, such that 

f{z) < y^, y~> 00. 

Note that it is relatively easy to generate Maass forms as eigenfunctions of the Laplace 
operator on a fundamental domain of the group F. 

Definition 2.17 A Maass form f is said to be a cusp form if the equality 

f f{z + b)db^O 
Jo 

holds for all z G H. 
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A Maass form / is call odd if / (— x + iy) = —f {x + iy) , and even if / {—x + iy) 
f{x + iy). 

One example of a Maass form is provided by the Eisenstein series: 



(m,n)5t(0,0) ' ' 



( 



= ^-^r(.)c(2.s) (I™ (7^))% 

7eroo\r 

where Too = {76 r|7oo = 00} . (9) 

(The factor 7r~*r (s) is introduced to get a simpler functional equation.) It can be checked 
that this series is an eigenfunction of the Laplace operator with eigenvalue s (1 — s) , and it 
turns out that it possess a nice functional equation. Eisenstein series can also be defined for 
different subgroups F and different cusps. 

From the spectral point of view, Eisenstein series correspond to the continuous part of the 
spectrum of the Laplacian. There are also other forms that correspond to the discrete part of 
the spectrum. 

By expanding a Maass form in Fourier series and taking the Fourier coefficients as the 
coefficients of a Dirichlet series, one can construct new L-functions that satisfy a functional 
equation. The construction is similar to the case of modular forms. 

Theorem 2.18 Let f be a Maass form with eigenvalue 1/4 + 1^^. Let e = or 1 depending on 
whether f is even or odd. Let 

A (s, /) = TT ^F ( \ r ( \ L (s, /) . 

Then A (s, /) is an entire function that satisfies 

A(s,/) = (-l)^A(l-s,/). 

Since the Laplace operator commutes with Hecke operators, therefore they can be simul- 
taneously diagonalized. 

For example, the Eisenstein series i? (z, i + it) is an eigenfunction of all Hecke operators 
with an eigenvalue 

„w=E(r- 

ad—n 

The Maass forms which are eigenfunctions of Hecke operators lead to L-functions with a 
Euler product property and in many respects the theory is parallel to the theory of L-functions 
which are obtained from modular newforms. 

In particular, for the L-series corresponding to the Eisenstein series LJ^ — iz, it -I- ir) , 
one can check that 

L{s,Er) =C{s + ir)C{s-ir). (10) 
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2.3.3. Statistical properties of zeros: Global scale 

■ 

where i]j is the Maass form with eigenvalue A and L is the corresponding L-function. 5* (i, '0) 
is related to the number of zeros of L in the critical strip in the same way as the usual S (t) 
function is related to the number of zeros of Riemann's zeta function. 

For this function S {t, tp) , Hejhal and Luo [16] proved analogues of Selberg's results for the 
function S (t) . 

For the Eisenstein series Er = E (^z , ^ + ir^ , formula (10) implies that 

S{t,Er)^S{r + t)-S{r~t), 

and the results about the statistical properties of S {t, Er) are reduced to the question about 
the statistical properties of S (t) . 

For the case of Maass cusp forms with eigenvalues from the discrete part of the spectrum, 
Hejhal and Luo proved the following result. 

Let Sj {t) S {t, tpj) where tpj has an eigenvalue A_,- = j + r^. Since ipj {z) is periodic along 
the real axis, it has a Fourier series, and it is known that it can be written as follows: 

{x + iy) = Y^ p, in) ^K,r, (27r \n\ y) e^""-, 

n#0 



where K^, [z) are the standard Bessel functions. Let Vj (n) = pj (n) / ^cosh irrj . The numbers 
i^j (1) will be used as weights in the limiting procedure. One knows that they are O {y/rj) and 

Since by the Weyl law, the number of rj below T is proportionate to T^, these weights can be 
thought as having bounded magnitude and not too sparse. 

Now we formulate one of the results in this paper (Theorem 3 in [16]). 

Theorem 2.19 (Hejhal-Luo) Let h > and t > be fixed. Then we have 

T^^ i2HT) ^^^^^^^ 2 (loglogT)"/^ 

where Cn are moments of the Gaussian distribution. 

The main new ingredient in the proof is the following lemma (Lemma 1 in [16]). 
Lemma 2.20 Let T > S, l3 > 0, ^ > e > 0, and let max(ni,n2) < T^"'^. Then 

oo „ -3/2 

J2e~PiT-r,)\^ (ni)i;~(;^= ^^T(5„,,„, +0^,, (d(ni)d(n2)(logT)2) , 
where d (n) denotes the number of divisors of n. 

The proof of this Lemma is based on the so-called Kuznetsov trace formula. 
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2.3.4. Local scale 

Rudnick and Sarnak [37] extended the results of Montgomery to zeta functions that arise from 
modular forms and from Maass forms. In fact, they work in greater generality and study the 
zeta functions that arise from the general automorphic cuspidal representations of GL,„ . The 
Hecke modular L-functions correspond to the case m — 2. Their main tool is the following 
explicit formula, which we formulate for the case of the Hecke L-functions. Let 



p\N ^^'^ p\N 



1 



c{p)p + p 



,2fe-l-s 



where 



Lp (s,/) = 



(1 - ai{p)p-'') (1 - a2(p)p-0' 

with the convention that for p\N one of (p) is zero. Let a (p'^) = ai (p)'' + a2 {p)^ , and 
define 6 (n) = A (n) a (71) . Then 



L' 



bin) 



n=l 



Theorem 2.21 (Rudnick and Sarnak) Let h e 

function, and let h(r) — J h (u) e"'"c?u. Then 



be a smooth compactly supported 



E'^W = ^ l°lh{r) (logg + E 



r /I 



r 2+^^ 



F' /I 



dr 



^/i(logn) + ^/i(-log7i) 



(11) 



where fij are some parameters that depend on the form /, and Q is the conductor of the form. 

By using this result and estimates on the size of coefficients b (n) , Rudnick and Sarnak 
proved a generalization of the Montgomery theorem. Their result is valid not only for the 
Riemann zeta function, but also for Dirichlet L-functions, for Hecke modular L-functions and 
for L-functions that correspond to automorphic cuspidal representations of GL^. We formulate 
it for Hecke modular L-functions. 

Consider the class of smooth test functions F{xi, . . . ,Xn) that satisfy the following condi- 
tions: 

TF 1. F{xi, . . . ,Xn) is symmetric. 

TF 2. F{x + i (1, . . . , 1)) = F (x) for aU t G M. 

TF 3. F (x) — > rapidly as |a:;| — > 00 in the hyperplane ^ Xj = 0. 

If Bn is a set of numbers xi, . . . , a; at, then the n-level correlation sum is defined by 



SCBf, 
\S\=n 
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Define the n-level correlation density by 

sin TTX 



Wn {Xl,.. .,Xn) = det {K {Xi - Xj)) , K (x) 



TTX 

Theorem 2.22 Assume the Riemann hypothesis for the zeros of L (s, /) . Let F {xi, . . . , xat) 
satisfy TF 1, 2, 3 and in addition assume that F [S^) is supported in < 1. Then, 

i?„ (Bjv, F) -> j F{x) Wr, {x) 5 j^^i±__±^^ _ _ 

as N ^ oo. 

Rudnick and Sarnak mention that the result can probably be proven for functions F with 
the Fourier transform supported in | < 2 by an improvement of their method, and 

conjecture that it holds without any assumption on the support of F (^) . 



2.4- Elliptic curve zeta functions 

The main source for this section is the book [30] by Milne. Consider an elliptic curve 

E : Y'^Z ^X^ + aXZ^ + bZ^, 

where a and b arc integer, and assume that |A| = 140"^ + 276^ | cannot be made smaller by a 
change of variable X ^ X/c'^, Y ^ Y/c^. This equation is called minimal. The equation 

E : Y^Z = X^ + aXZ^ + 

with a and b the images of a and b in Fp (the finite field with p elements) is called the reduction 
of E modulo p. (It is assumed here that p 7^ 2, 3. In the case when p is 2 or 3, a somewhat 
different notion of the minimal equation is needed.) Let Np is the number of solutions of this 
equation in Fp. 

There are four possible cases: 

(a) Good reduction. E is an elliptic curve. (That is, the determinant does not vanish and 
therefore the curve is smooth.) This happens if p 7^ 2 and p does not divide A. In this case we 
define the number of points of E by Np. 

(b) Cuspidal, or additive, reduction. This is the case in which the reduced curve has 
a cusp. For p 7^ 2, 3, this case occurs exactly when p\Aa^ + 27&^ and p\ — 2ab. 

(c) Nodal, or multiplicative, reduction. The reduced curve has a node. For p 2,3, 
it occurs exactly when p\Aa'^ + 27b^, p | —2ab. 

(cl) Split case. The tangents at the node are rational over Fp. This happens when 
—2ab is a square in Fp. 

(c2) Non-split case. The tangents at the node are not rational over Fp.This occurs 
when — 2a& is not a square in Fp. 

The names additive and multiplicative refer to the group of points on the reduced curve, 
which in these cases is isomorphic either to (Fp, +) , or (F*, x) . 

We define the L function associated with the elliptic curve E as follows. 
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Here, the local factors Lp (T) are defined as follows: 



apT + pT , if p is good, with ap = p + 1 — Np, 

1 — T if i? has split multiplicative reduction, 

1 + T if has non-split multiplicative reduction, 

1 if £' has additive reduction. 



Let S be the (finite) set of primes with bad reduction. Then we can also write 

^^''^^ ' ^ n mf^ n i+(Np-p~i) + 



pes -F V/' ; 



n Lp ip-^) n (1 _ app-s) (1 _ ^pp-s) ■ 



The Hasse-Weil conjecture says that L (s, E) can be analytically continued to a meromor- 
phic function on the whole of C and satisfies a functional equation. The recent work by Wiles 
and others confirmed this conjecture by showing that all elliptic curves are "modular" , in 
particular, their L-functions arise from modular forms. To a certain extent, this result reduces 
the study of the elliptic L-functions to the study of the Hecke modular L-functions. 

It is known that the numbers ap do not exceed in absolute value. For a fixed elliptic 
curve and different primes p, these numbers are believed to be distributed on the interval 
[— 2y^, 2^y/p] according to the semicircle distribution but this is not proven. In fact, this 
conjecture is related to the Birch-Swinnerton conjecture that states that 

L{s,E) - C(s- 1)'' as s ^ 1, 

where r is the rank of the group of rational points on E, and C is a certain predicted constant. 
A stronger statement, which can still be true, is that 



Pis pts 

p<x 

It is possible to construct zeta functions for other nonsingular projective varieties and the 
conjecture by Hasse and Weil states that these zeta functions satisfy a functional equation 
(and the Riemann hypothesis). However, apparently not much is known beyond the cases of 
projective spaces and elliptic curves. 

3. Selberg's zeta functions for compact and non-compact manifolds 

The main source for this section is Hejhal's book [15]. 

3. 1 . Selberg 's zeta function and trace formula 

Let M be a compact Riemann surface of genus g > 2. Then, M can be identified with a 
quotient space r\H, where H is the upper half-plane and F is a discrete subgroup of SL2 (R). 
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We assume that H has the Poincare metric \dz \ /y with the area element dxdy/y^, and therefore 
it has the constant negative curvature. This metric is naturally projected on the surface M. 

This is not the most general situation of interest since most of the quotient spaces r\H 
occurring in arithmetic applications have cusps and therefore are non-compact. However, the 
theory is most clear and transparent for the compact surfaces. 

The Laplace operator on M can be defined by the following formula. 

-A : u -y"^ {uxx + Uyy) . 

It can be shown that this operator has a discrete set of non-positive eigenvalues: 

= Ao < Ai < A2 < . . . , 

and the only point of accumulation of these eigenvalues is 00 . 
Let us define 

An — -J , if A„ > -J , 



A„ + i, if A„ < 1/4, 

so that A„ = ^ + r^. 

Also let m = max{fc : A^ < 1/4} . 

Let [M) be the set of all closed geodesies on Af , and let V (M) be the subset of all 
prime closed geodesies (that is, the closed geodesies that cannot be represented as a non- 
trivial multiple of another closed geodesic). It is known that it is a countable set, which we 
can order by the lengths of its elements. Closed geodesies correspond to hyperbolic elements 
of the group F (that is, the elements of F with the trace outside of [—2, 2]) up to conjugacy of 
these elements. If P G F corresponds to a geodesic 7, then 7 is prime if and only if there is no 
Pq € F such that P = Pq' for an integer fc > 1. 

If I (7) denotes the length of the geodesic 7, then we set 

I7I :=e'(''). 

If P is the hyperbolic element of F that corresponds to 7, then 

|7|V2 + |^|-i/2^|TrP|. 

We will also write I7I = N [P] (meaning norm of P). 

The Selberg trace formula relates sums over eigenvalues A^ to sums over hyperbolic el- 
ements (geodesies) [P]. Let h{u) be a function which (i) is analytic in the strip |Imu| < 
1/2 + 5, (ii) is even: h{u) = h{—u), and (iii) declines sufficiently fast in the strip: |ft.(u)| = 



0((1 



2 S\ 

\Rcu\y^'°]. This ensures absolute convergence in the following formulas. 



Let h{t)^ ^Jh (u) e-**"dM. 
Then 



E 

n=0 



, h{rn) = ^-^—^ / r/i (r) tanh (7rr) dr (12) 
- 27r Jr 
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where the sum is over all distinct conjugacy classes of hyperbolic elements [T] and [Tq] is the 
primitive element for T,T ^ Tq. 

It is instructive to compare this formula with formula (11). Since it resembles the explicit 
formulas from number theory, it is natural to seek for a corresponding zeta functions. 

We define the Selberg's zeta function as follows: 



It turns out that the Selberg's zeta function is closely related to the eigenvalues of the 
Laplace operator on M. 

Theorem 3.1 (Hejhal-Selberg) (a) Z (s) is actually an entire function; 
(h) Let P he a real number > 2. For all s, the following identity holds: 



where /i (F) is the area of the fundamental region of the group T. 

(c) Z (s) has "trivial" zeros s = — fc, fc > 1, with multiplicity (2g — 2) (2fc + 1) ; 

(d) s = is a zero of multiplicity 2g — 1; 

(e) s = 1 is a zero of multiplicity 1; 

(f) the nontrivial zeros of Z (s) are located at ^ zt zr„. 

Since all but a finite number of eigenvalues are greater than 1/4 hence all but a finite 
number of r„ is real and therefore the claim (f ) implies that all but a finite number of zeros 
of Z (s) are located on the line Imz = 1/2. 

The formula in claim (b) of this theorem follows from Selberg's trace formula and it can 
be thought as a functional equation for the logarithmic derivative of Z (s). In particular, it 
implies the functional equation for the zeta functions itself. 

Theorem 3.2 (Hejhal-Selberg) Selberg's zeta function satisfies the following functional equa- 



oo 




(13) 



jeV(M) k=0 




tion: 




3.2. Statistics of zeros 



The number of zeta zeros in a long interval has the following asymptotics: 



N[k -.Q <rk <T 



1^{F) 



rp2 



+ S{T) + E (T) , 



where 
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and 



E{T)^0 (1) = 2c / t [tanh (7ri) - l]dt - {m + I) . 







In other words, the number of zeros in a unit interval is ^ cT. In comparison, for Riemann's 
zeta we have ^ c log T zeros in the unit interval. 
It is known that 



S{T) = 



T 



logT, 



and S (T) = fl± 



1 \ 1/2' 

log 1 ^ ' 



log log T 



(Recall that the notation / (x) = fi+ [g(x)] means that limsup > 0, and / (x) = 

r2_ \g{x)] means that liminf (^^^^ < 0-) 

It was found that it is difficult to generalize the results concerning the moments of the 
S [x) function to the Selberg's zeta. Since these results are essential for the study of statistical 
properties of zeta zeros, there is a stumbling block here. 

Selberg managed to generalize some of these statistical results for a particular choice of the 
group r. 

Let p > 3 be a prime and A be a quadratic non-residue modulo p. Define 




1=1 {A,p) = < rr 2/0,2/1,2/2,2/3 are mteger. 

[ V 2/2VP- 2/3V-4p 2/0 

and call it a quaternion group. 

It turns out that F {A,p) has no parabolic elements, (that is, it has no group elements with 
trace equal to ±2, except ±7). Moreover, it corresponds to a compact Riemann surface. 

If p = 1 (mod4) , then F {A^p) has no elliptic elements. 

Let S (t) = S+ (t) - S- (t) , where S+ (t) = max{0, S (t)} and (t) = max{0, -5 (t)} . 
Then the following theorem holds. 

Theorem 3.3 (Hejhal-Selberg) Let F = F {A,p) with p=l (mod4) . Then (for large T): 

1 r'^^ T 
S+ {tfdt > ci- 



tJt ^ '^""■'^logrf 

where ci is a positive constant that depends only on F. A similar inequality holds for (t) . 

In order to appreciate this result note that it implies that the average deviation of S (T) 
from its mean (if it exists) is of the order larger than Vt/ (logT) which should be compared 
with the number of zeros in the interval [0,T] , that is, cT^. To put it in prospective note 
that the average deviation of the zeros of S (T) for Riemann's zeta function is of the order 

1 /2 

(log log T) which should be compared with cT log T, the number of Riemann's zeros in [0, T] . 
These situations appear to be quite different. 

Moreover, recently there was some numeric and theoretical work on the eigenvalues of 
the Laplace operator on manifolds T\II for arithmetic groups F. First, numeric and heuristic 
analysis showed that the spacings between eigenvalues resemble spacings between points from a 
Poisson point process rather than spacings between eigenvalues of a random matrix ensemble 
(see Bogomolny et al [5] and references wherein). Next some rigorous explanations of this 
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finding liave been given. They are based on exceptionally large multiplicities of closed geodesies 
with the same length. See Luo and Sarnak ([27] and [28]) and Bogomolny et al [6]. 

There is also some work on correlations of closed geodesies - see Pollicott and Sharp [34] . 

3.3. Comparison with the circle problem 

If we consider the compact Riemann surface of genus 1, then we are led to similar questions. 
Such a surface can be represented as quotient space A\C, where A is a lattice. Consider, for 
concreteness, A = . Then the eigenvalues of the Laplace operator are 47r^ (m^ + n^) , 
where m and n are integer, and the number of the eigenvalues below t equals the number of 
integer points in the circle t/n. Let 

r {n) = N {{a,b) e Z X Z : a'^ + b'^ = n} , 

and 

A (x) ~ r (n) = TTX + R (x) . 

Q<n<x 

The function A (x) can be thought as the counting function both for eigenvalues of Laplace 
operator and for closed geodesies of bounded length. 

Then by using the Poisson summation formula it is possible to derive the following result: 

POO 

y2 r (??.) / (n) = TT V] r (??.) / / (x) Jq {2T:y/nx) dx. 

Jo 

Given that it is possible to apply this to the step function / (x) , one can obtain the explicit 
formula: 

oo , s 

R(x) = y/^y^—^Ji (27rv^) . 

n— 1 ^ 

This leads to various estimates on R (x) , in particular it is known that 

R = (x^^^) and R = n± (x^^'^) , 



and that ^ 

- [ R {tf dt = + o [(log xf 



X 







This suggest that the "standard deviation" of R{t) is x^/^. Again this does not resemble 
the situation with Riemann's zeta zeros. 

Some more details about this problem and references to the early papers can be foimd in 
[24]. More recent research include [14] and [4]. 

4. Zeta functions of dynamical systems 

The main sources for this section are reviews by Ruelle ([38] and [39]) and Pollicott ([36] and 
[35]). 
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4-1. Zetas for maps 

Let / be a map of a set M to itself, let the periodic orbits of / be denoted by P, and let \P\ 
denote the period of orbit P. Then, we can define the zcta of / by the following formula: 

c(z) = n(i 
p 

= exp^— |Fix,n, 

m— 1 

where |Fix /™| denote the number of fixed points of /™. 
Permutations 

Let M be a finite set, and let / be given by a permutation matrix A. Then the number of 
fixed points of /™ is given by TtA™. Hence, we have 

C{z) = expTr^^"^) 



m 

rn— 1 

= exp(-Trlog(l - z^)) 
= l/det(l-z^), 

which is closely related to the characteristic polynomial of matrix A. In particular, all poles 
of the zeta are on the unit circle. 



4- 1.2. Diffeomorphisms 

Let M be a torus R^/Z^ and let / be induced by a linear transformation A G SL2 (Z) . Assume 
that the eigenvalues of A are positive and not on the unit circle: Ai > 1 > A2 > 0. Then the 
number of the fixed points of /™ is the number of solutions of the equation / — A"^ = 
(modZ^) , which can be computed as |det (/ — A"^)\ = (A™ — 1) (1 — A™) . Hence, we have 

00 

C{z) = expTr V -[(zAi)™ + (zA2r-z'"-(zAiA2)n 
^-^ m 

m—l 

(l-zAi)(l-zA2) 
(l-z)(l-zAiA2) 
(l-zAi)(l-zA2) 
{l-zf 

This example can be generalized to linear transformations A e GLd (Z) acting on a d- 
dimensional torus or more generally to discrete groups acting on compact quotients of nilpotent 
Lie groups (nil-manifods). 

If we wish to generalize this example to diffeomorphisms of smooth orientablc manifolds, 
it is natural to count the fixed points by taking into account the degree of map / at the fixed 
points. Recall that the degree equals to 1 if the map is non-singular and preserves orientation 
at the fixed point, and to —1 if it reverses the orientation. It can also be generalized to the 
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situation when the map is singular at the fixed point. If d{x, f) denotes tlie degree of / at a 
fixed point x, then we define the Lefschetz zeta function as 



a(z) = exp5]— d{x,r). 



m 

ni=l a:eFix(/'") 

In this case one can use the Lefschetz fixed point formula that says: 

dim M 

a;GFix(/'") i=0 

where Hi is the i-th homology group of the compact manifold AI with real coefficients, and 
is the map induced by on Hi. 
In particular, if Xij are eigenvalues of then we get 

dim M dim Hi 



m— 1 i—0 J 

dim M I dim Hi 

= n n (1^^.) 




t=0 \ j=l 

which is a rational function. 

If / is a diffeomorphism of a compact manifold M, then the original dynamical zeta (counted 
without taking into account the degree of Z™) is often called the Artin-Mazur zeta function 
(see Artin-Mazur [2]). If the diffeomorphism satisfies some additional conditions (hyperbolicity 
or Axiom A), then it is known that this function is rational. (This was conjectured by Smale 
[45], and proved by Guckenheimer [11] and Manning [29].) 



4.1.3. Subshifts 

Suppose next that A is an N-hy-N matrix of zeros and ones, and that the set M consists 
of doubly infinite sequences {xi} of symbols 1, . . . , A^, which satisfy the following criterion. A 
sequence {xi} belongs to M if and only if A^iXi^i = 1 for every i. In other words, the symbol 
Xi determines which of the other symbols are possible candidates for Xi+i. The map / is simply 
a shift on this set M: {xi} — {xi+i} . In this case, the number of fixed points of J™ is simply 
Tr (A™) , and we have 

e(2) = l/det(l-zA), (14) 

similar to Example 1. Note that if A is symmetric then the dynamics of the system is in a 
certain sense reversible. On the other hand the poles of the zeta in this case are all real. So 
this example illustrate a connection between the distribution of zeros and the reversibility of 
a system. 



4-. 1.4- Ihara-Selberg zeta function 



A basic reference for this section is a book by Terras [46]. 
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Let G be a finite grapli. Orient its edges arbitrarily. Let the 2\E\ oriented edges be denoted 
, . . . ,62,1 = e„ . A path is a sequence of oriented edges sucli that 
the end of one edge equals to the beginning of the next edge. A path is closed if the end of 
the last edge corresponds to the beginning of the first edge. A closed path (ei, 62, . . . , es) is 
equivalent to the path (e2, ...,6^,61). A closed path P is primitive if P 7^ D"^ for m > 2 and 
any other path D. A path is non-backtracking if e^+i 7^ e^^ for any i. A closed path is tailless 
if fls =/= a^^ . 

If P = (ei, 62, . . . , 65) , then I (P) = s is the length of the path. 

A prime [P] is an equivalence class of primitive closed non-backtracking tailless paths. (It 
is an analogue of a geodesic on a manifold.) 
The Ihara zeta function is defined as follows: 

[P] 

It has a representation as a determinant which was first discovered by Ihara for regular 
graphs ([19]) and then by Bass [3] and Hashimoto [12] for arbitrary finite graphs: 

Cg (u)-' = (1 - u') dct (/ - Agu + Qgu') , 

where Ag is the adjacency matrix of G, and Qg is the diagonal matrix whose j-th diagonal 
entry is (— 1 + degree of j-th vertex). 

The Ihara zeta function can be though as a particular case of a dynamical zeta function 
for a subshift associated with graph G. The alphabet consists of 2 \E\ oriented edges, and the 
subshift matrix equals to the edge adjacency matrix Wg which is defined as follows. 

Definition 4.1 The edge adjacency matrix Wg is the 2 \E\-by-2 \E\ such that i,j entry equals 
1 if the terminal vertex of edge i equals the initial vertex of edge j and edge j is not the inverse 
of edge i. 

In particular, from (14) we have another detcrminantal formula: 

(g(u)-i ^ dct {I -uWg). 

If the graph is g + 1 regular, that is, if every vertex has degree q + 1, then there is a pole 
of the zeta at u = 1/q. Moreover, every non-real pole of the zeta lies on the circle with radius 
\j ^fq. The Riemann hypothesis for regular graphs says that all non-trivial poles are on this 
circle. This is not always true, and a useful criterion is as follows. Let the second largest in 
magnitude eigenvalue of Ag be denoted Ai. Then the Riemann hypothesis holds if and only 
if I All < 2y/^. The graphs that satisfy this condition are often called Ramanujan following 
a paper by Lubotsky, Phillips, and Sarnak [26], which constructed an infinite family of such 
graphs. 

Note that 2^ is critical in the sense that for every e > 0, there exists a constant C such 
that for all (g + l)-regular graphs with n vertices there is at least Cn eigenvalues in the interval 
[2y^ — e, g + 1] . On the other hand, a random regular graph is approximately Ramanujan 
with high probability. That is, for arbitrary e > 0, the probability that |Ai| > 2y^-|-e becomes 
arbitrarily small as the size of the graph grows. This fact is known as the Alon conjecture and 
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it was recently proved by J. Friedman in Friedman (2008). However, the distribution properties 
of the largest eigenvalue are still unknown. 

In his Ph.D. thesis [33], Derek Newland studied spacings of eigenvalues of random regular 
graphs and spacings of the Ihara zeta zeros and found numerically that they resemble spacings 
in the Gaussian Orthogonal Ensemble. See also an earlier paper by Jacobson, Miller, Rivin 
and Rudnick [22]. 

For Ihara'z zeta function, there is an analog of Selberg's trace formula which we formulate 
for the case of a g + 1 regular graph. 

Let Nm be the number of all non-backtracking tailless closed paths of length m. (Note that 
this is the number of all closed paths, not the number of the equivalency classes.) In particular. 



du 



This shows that 



[P] 

E 

p 



Cg (u) = cxp Nm — . 

— ' m 



m>l 

Theorem 4.2 (Terras) Suppose < a < 1/q. Assume that h{u) is meromorphic in the 
plane and holomorphic for |u| > a — £, e > Assume that h (u) = O (^\u\ , a > 0. Let 

h (n) :— 2^ J\u\=a ('^) '^^ '^"■'^ assume that h (n) decays rapidly enough. Then, 

^ ph (p) = Nnh (n) , 

p n>l 

where the sum on the left is over the poles of C,q {u) . 
Some applications of this result can be found in . . . 

4.1.5. Frobenius maps 

Another important example arises if M is the set of solutions of a system of algebraic equations 
over the algebraic closure of ¥q (i.e., the finite field with q = p"^ elements), and / is the 
Frobenius map Frob: a; — > x'. These functions are often called the Weil zeta functions although 
it should be noted that the first examples were originally considered by E. Artin. 

Weil's zeta functions directly connected to the number-theoretical zeta-functions. Consider 
the affinc curve given by the equation f {X,Y) = over the finite field F^. Let p denote a 
prime ideal of the field F, [X, Y] // {X, Y) and let the order of F, [X, Y] /p be denoted by Np. 
Then, by analogy with Ricmann's zeta function we can define 



It turns out that this definition is in agreement with another definition of the Weil zeta function: 

— S7TL 

C (s) = exp Y Nr. ^ 



m 
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where Nm is the number of fixed points of Frolf^ , or to put it in other words, the number of 
solutions of the equation / (X, F) = in the finite field F^m . 

Here is an example. Consider the affine line. Then = g™, and 

C{z) = expVA'"^ — 
^-^ m 

m— 1 

= exp [- log (1 - qz)] 
_ 1 
1 — qz 

In another example, the curve is elliptic. Note that the Frobenius map is a linear operator over 
Fg and there is some resemblance to Example 2 about linear maps on a torus. We can guess 
then that 

1 - zTr [Frob] + det [Frob] 
^ " (1 - z) (1 -z det [Fro5]) ' 
provided Tr [Frob] and det [Frob] are given some appropriate meaning. If we believe det [Frob] = 
q, then it A^i = det (/ - Frob) = det(_Fro6) - Tr {Frob) + 1, and so Tr [Frob] = q + l - Ni. 
Hence, our guess would be 

l-z(g + l-iVi) + z^g 
(l-z)(l-zg) ' 

and it turns out that this guess is correct. Of course. Ad is in fact infinite-dimensional over 
Fq, so we do not know what is the meaning of the determinant and the trace. However, one of 
the proofs of the rationality of the zeta function uses a two-dimensional module over Q; (the 
field of ?-adic numbers) and the arguments resemble our heuristic. One can consult the book 
by Silverman [44] to see more. A more standard proof is based on the Riemann-Roch theorem 
and counting divisors. It can be found in [30]. 

The Riemann hypothesis in this example is equivalent to the statement that 

[N^-q-\[<2^, 

since this implies that the roots of the polynomial in the numerator are on the circle with 
radius 

Artin [1] proved the zetas of the curves can be explicitly computed as rational functions of 
z = p^". He checked numerically that in many cases the zeta function satisfies an appropriate 
variant of the Riemann hypothesis, and conjectured that this should be true for all curves. 
Later this was proved by Hasse [13] and Weil ([50] and [49]). 

Weil also conjectured in [48] that the zeta function of every algebraic variety is rational, 
that it has a functional equation, and that it satisfies the Riemann hypothesis. The proof of 
this general conjecture led to an introduction of many new ideas in algebraic geometry by 
Dwork, Grothendick, and Deligne. 

Katz and Sarnak [23] have studied the distribution of the zeros of the Weil zeta functions 
when the genus of the corresponding curve grows to infinity. They found that for "most" of the 
curves the local statistical behavior of these zeros approaches the behavior of the eigenvalues 
from the random matrix ensembles. It is interesting that there are no explicit constructions of 
curve families for which it is known to be true. It is only known that it must be true for most 
of the curves. 
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4- 1.6. Maps of an interval 

For yet another example consider the map 2; — > 1 — ^x'^ of the interval [—1, 1] to itself. For a 
speeial value of « 1.401155... (the Feingenbaum value), this map has one periodic orbit of 
period 2" for every integer n > 0. Therefore, 

n=0 Ti=0 

This C satisfies the functional equation C (z^) = (1 — z) ^ (z) . More generally, the piecewise 
monotone maps of the interval [—1,1] to itself correspond to the Milnor- Thurston zeta func- 
tions. See [31]. As far as I know there was no systematic study of statistical properties of their 
poles and zeros. 

4-2. Zetas for flows 

If / is a flow on M, that is, a map M x IR+ AI, then we can define the zeta function of this 
flow as 

c(.)=n(i-e-^'^"0"' 

UJ 

where w denotes a periodic orbit of /, and / (w) is its length. It is a more general case than 
the ease of maps since there is a construction ("suspension") that allows us to realize maps 
as flows but not vice versa. Unsurprisingly, it turns out that zeta functions for flows are more 
difficult to investigate than zeta functions for maps. 

If we imagine that prime numbers correspond to periodic orbits of a flow and that the 
length of the orbit p is given by logp, then the zeta function of the flow will coincide with 
Riemann's zeta function. 

One particularly important example of a flow is the geodesic flow on a smooth manifold 
M. In the case when M has a constant negative curvature, the corresponding dynamical zeta 
function is closely related to the Selberg zeta function. Namely, 

C is) = ^yT;^ ^^^^ Zis)^f{Cis + n)-\ 

^ ' n=0 

Another important example is the geodesic flow for billiards on polygons. 

I am not aware about the functional equation for zeta functions that comes from general 
geodesic flows (other than flows on manifolds of a constant negative curvature). Also, it 
appears that not much is known about the statistical properties of the distribution of zeros 
and poles of these zeta functions. 

The main benefit of these functions is that the distribution of closed geodesies depends on 
the maximal real part of their poles, and this quantity can be studied by methods that do not 
involve the functional equation or trace formula. Hence, one can study geodesies on spaces of 
variable curvature. The idea of this approach is that geodesic flows can be written as shifts 
of sequences of symbols and poles of zeta functions are related to eigenvalues of operators 
induced by these shifts. (This is the method of symbolic dynamics). 
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5. Conclusion 

We considered tlie statistical properties of the zeta functions. In summary, the zeros of number- 
theoretical zeta functions satisfies many properties which are true for eigenvalues of random 
matrices. The main outstanding problem is the Montgomery conjecture that says that the 
correlations are the same as for random matrices even if the test functions are not smooth. 

The statistical properties of the Selbcrg zeta zeros are different from the properties of the 
eigenvalues of random matrices. It is conjectured that the local properties coincide with the 
properties of the random independent points However, the exact description of these properties 
is not known. 

The zeros of dynamical zeta functions are often related to eigenvalues of certain matrices, 
and the study of their properties depends on the study of the matrix properties. Two interesting 
cases are the case of zeta functions of algebraic curves and the case of Ihara's zeta functions 
of graphs. 
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